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Within the frame-work of semiclassical theory two-level approximation in atomic system has 
been considered. Model proposed by M.D. Crisp and E.T. Jaynes has been modified. It has been 
shown that the time-dependent frequency shift depends on the higher multipole moments, retained 
in the Taylor expansion of electromagnetic field. 
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Issues, related to the problem of interaction of photon and micro-particle, in their full length are beyond the scope 
O ' of quantum mechanics. They cannot be considered without invoking additional principles concerning the laws of 
f \ \ occurrence and disappearance of electromagnetic field. According to quantum mechanics atom should remain in 
excited state for long in absence of external field, whereas experiment shows that atom transforms into normal state 
emitting photon. This contradiction can be explained if we take into account the fact that the moving electron creates 
electromagnetic field which acts on the electron. Several authors tried to consider this reverse action of field on 
electron several ways. One of these methods was proposed by Jaynes and Cummings jl] in 1963 that was further 
k*" . developed by Jaynes and Co and many others ■ 

In classical electrodynamics the radiative process are calculated from self-energy of the electron in external fields. 
In contrast, in quantum electrodynamics, the self-energy is first thrown away and one begins with bare particles; then 
the self-energy is put back in photon by photon, hence the use of perturbation theory. Recently, Barut and coauthors 
developed a quantum electrodynamics based on self-energy ||||]. 

■ Authors of the papers mentioned previously mainly confined their study within electric dipole moment. Here 
' we make an attempt to enlarge this study taking into account the moments of higher order, particularly electric 
, quadrupole moment. 

1 Let us consider a nonrelativistic, spinless particle in external magnetic field. It can be described by the Hamiltonian 



(N 



^1" 
O 



ij = ^r p _£ A l 2 -- (i) 

C3 2m l - J 

'■ Varying this Hamiltonian with respect to A and using the continuity equation + divj = one finds 

j = !£^{^v^* - - — A1"5* (2) 

rN , 2m mc 

p = evTvP (3) 

Taking the field to be weak one we further neglect the diamagnetic term in the Hamiltonian and current density. 
Now, any state of atomic system may be expressed as 



tf(r,*)=53a a (t)lMr) (4) 



where 4>( v ) is the eigen functions of Hq = — (h 2 /2m)\7 2 — e 2 /r, i.e., 

HoV'a(r) = E a ip a (r) 

Putting (4) into (3) we obtain 
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where 



1 



Bijan Saha 



P0a(t) = a a (t)a*p(t) = p a p(t)* 
is the f3a density matrix element of the atom in the Schroediger picture that evolves according to 
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Since the magnetic field obeys the Maxwell equations, for A in Coulomb gauge (divA = 0) we can write 
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„. 1 d 2 A 4w. ± 
V 2 A- -— = j x 



c 2 dt 2 

Here j is the transverse current density and defines as 

= — V x V x 
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j(*,x') 



d 3 x' 



Further we denote j = j. The solution to the Maxwell equation can be written as 

A(x,i) 



1 /-j fr'.t-lx-x'l/c) ^ 



Taylor expanding the expression for j one gets 

j(x',*) j3 , i rdj(x',t)„ , , i /-a 2 j(x',i) 



A(x, t) 



r <fx' - - 
|x — X' c 



dt 2c 2 



x-x'|d 3 x' 



dt 2 



Further expanding |x — x'| for x' « x one finds 



A(x,t) 



•?( x '>*)_ d 3 x / _ I /" 9 i( x '>*) j3„/ , a; fd 2 j(x',t)^ r/ x 



dt 



<% 2 



2c 3 r 



(9) 



(10) 



(11) 



Putting \& = ^ Q a a (t)ip a (x., t), where ^: Hoip a = E a ip a into the equation above and retaining the electric dipole and 
quadrupole moments we find 



A(x,i) «5^/>o0(t) 

a/3 





f dk J 


-2ir 2 mc J 





dk I dn(P\e- ik - x V\a) ± e lk -' x 
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where the transition frequencies and the electric dipole and quadrupole moments are defined, respectively, as 

tt aP = {E a - Ef})/n, (13a) 
T) a p — J ^p a exippdx 7 or in components -D^g = J ^p a ex l ippdx (13b) 

Q% ] = J VaerVppdx, r« = i(.xV - V^'), r = |x|, 

Here Ao is an externally applied field. Putting the expression for A into (R), for density matrix we find 

Pa/3 = —i&a/3PaP ~ «^(F QK - V K/ g) p KK p a {3 

K 

- ^[^(^(W + - + Bp K ) + {C aK + Cp K )]p KK p a p 

K 

where we define 



(14) 
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r «/3 = - 2 ^2 C 2 / / dn(a\e* k - x ' |/3)j.(/3| e -*» ± = (15a) 
o 

A a p = ^(Do^D^/fic 3 )!!^ = — Ap a , Einstein coefficient (15b) 

B a p = (DapAap/hc^ttlp =, -Bp a , A af3 = J rJ a pW, (15c) 

/k 
^-J Q/3 (x)dx. (15d) 

Here we denote J a p = (efr/2mi) [ippVipa — 4>p V^ a ] ■ 

The equation (jlj) can be written in the following way where the repeating index denotes summation 

PaP = —iQaPPjrMapyr — i(T aK — T K p)p KK p 1T M a p 1T 

r -(A aK +Ap K ) - (B aK + Bp K ) + (C aK + Cp K )]p KK p 1T M a p 1T (16) 
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[Sl aK D$p KT - Cl K pD£) p 1K ]M a p 1 r, 



where M a p 1T — S a ^5g T . 

As one sees from (UJ) or (16), the off-diagonal density matrix elements oscillate at frequencies f!, a p + S£l a p(t), where 
the time-dependent frequency-shift is 

5Q a p(t) = - ^(r aK ~ r K p)p KK (t) (17) 

K 

Now the expectation of dipole moment of the atom 

<p>= J **(x,i)ex*(x, t)dx 

in account of (^) can be written as 

< p >= ^D Q(3 p ( g Q (i). 

a/3 

Thus we see that the off-diagonal matrix elements are directly connected with the expectation of dipole moment. 

In what follows we take into account only two of these levels. We choose the zero from which we measure the 
energies to be midway between the two active levels, so that 

E 2 = -E x (18) 

The equation (|l6| ) can then be written as 

p n = ~2qpnp 2 2 (19a) 
p22 = 2qpnP22 (19b) 

Pl2 = — * [Ol2 + Tiipn - T 22 P22 - r i2 (pll - P22)]P12 + q(pil ~ P22) Pl2 (19c) 

P21 = — * [0 2 i - Tnpii + T 22 p22 + r 12 (pn - p 2 2)]p2i + q{pn - P22)P2i (19d) 

where we denote 2q = A12 — 2B\2 + 2Cp. 
Let us now rewrite p a p in the form [|lO|]ll| l 

Pc t p = l(S a p + P j a 1 a0 ) (20) 

where a 3 are the Pauli matrices and P = (P x ,P y ,P z ) is a unit vector of three-dimensional Poincare representation. 
From (Sol) follow: 
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Pll = \(l + Pz), PX2 = \{Px-iPy), 



P22 = \{l~P z ), p 2l = \(P x +iP y ) 



or equivalently, 

Pll+P22 = l, Pit - P22 = Pz, Pl2+p21=Px, Pl2 - P21 

In account of (|2l]) and ( p2| ) from (|l^) we find the following system of equations 



-iP» 



P x = qP z P x + (0l2 + r + AP z )P y 

Py = qP Z Py - (fiia + T + AP^P* 

P = g(i* - I) 

where we denote r = (rn — T22)/2 and A = (T22 +rn)/2 — T^. The solutions to the system of equations 

P x = cos [SI12 — to) + T (t — to) + (A/g)lncoshg(t — to)] sechg(t — to) 
P y = sin [fii 2 (t — to) + T (t — to) + (A/g)lncoshg(i — to)] sechq(t — to) 
P z = — tanhg(t — to) 

Rewriting ( pij ) in terms of p we find 

p u = 1/ exp [2q(t - to)] + 1 
paa = l/[exp[-2g(t-to)] + l 



Pl2 
P 21 



exp ( — i[fli2(t - t ) + r(t - to) + (\/q)bx cosh q(t - t )]j sechq(t - to), 
exp fi[Oi2(t — t ) + r(t — t ) + (A/^lncoshq^t — t Q )]j sechg(t — t ). 



For the expectation value of the energy in account of (|l^) we find 

< Ho > = E lPll (t) + E 2 p 22 (t) = 
= -^0 2 itanh[g(t - to)]. 



-0 2 l(/?22 - Pll) 



where as, for the expectation of the dipolc moment we obtain 

< p > = D 2 i(pi2 + P21) = ~D 21 P X 

= D 2 i sechg(t — t ) cos \p, 2 \t + $(t)], 

where we define 

0(t) = i? - rt - (X/q)lncoshq{t - to), i?o = [(Tn - r 22 )/2 - Q 21 ]t 
and corresponds to a time-dependent frequency shift 

Sfl 2 x(t) = M/dt — — t — Atanhg(t - t ) 
Comparing ( |29| ) with those obtained in S one finds the additional frequency shift as 

tanh[(C 2 i - S 2 i)(t - t )]sech 2 [^2i(t - t )/2] 



A(Sfl 21 {t)) = A 



1 + tanh[^2i(t - t )/2]tanh[(C 2 i - B 21 )(t - t )] 



(21) 



(22) 

(23a) 
(23b) 
(23c) 

|) read 

(24a) 
(24b) 
(24c) 

(25a) 
(25b) 
(25c) 
(25d) 



Thus we see that beside Einstein A coefficient, higher multipole moments, in particular quadrupole one, 
to the spontaneous decay of the atom from an exited state. 
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